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Using data taken with the CLEO II detector at the Cornell Electron Storage 
Ring, we have determined the ratio of branching fractions: Ry = T(T(1S) — > 
7flfflf)/r(T(lS) -> ggg) = (2.75 ± 0.04(stat.) ± 0.15(syst.))%. From this ratio, 
we have determined the QCD scale parameter (defined in the modified 
minimal subtraction scheme) to be A-^r = 233 ± 11 ± 59 MeV, from which 
we determine a value for the strong coupling constant a s (M T ( 1S )) = 0.163 ± 
0.002 ± 0.014, or a s {M z ) = 0.110 ± 0.001 ± 0.007. 

PACS numbers: 13.40.Dk, 13.25.+m, 14.40.Jz 
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Typeset using REVTgX 



I. INTRODUCTION 



The three primary decay modes of the T(1S) are three gluons {ggg), a virtual photon 

999 ^ a s ' 



(7*), or two gluons plus a photon (jgg). We expect these decay widths to vary as: Y ggg oc a 3 



I\. oc ac m , and Y~ nn g, respectively. From the ratio of decay rates: 

p — ^199 „ ^T(lS)^7gg a em 

K y = - — oc — oc — — , (1) 

J- ggg ^X(lS)-^ggg «s 

one can determine a value for the strong coupling constant, a s . In this analysis we determined 
this ratio by measuring the number of direct photon and three gluon events from a sample 
of T(1S) data. 



II. DETECTOR, DATA SAMPLE, AND EVENT SELECTION 

The CLEO II detector is a general purpose solenoidal magnet spectrometer and calorime- 
ter. Elements of the detector, as well as performance characteristics, are described in detail 
elsewhere 0. For photons in the central "barrel" region of the Csl electromagnetic calorime- 
ter the energy resolution is given by 

where E is the shower energy in GeV. The tracking system, time of flight counters, and 
calorimeter are all contained within a 1.5 Tesla superconducting coil. 

The data used in this analysis were collected on the T(1S) resonance at a center-of- 
mass energy E cm = 9.46 GeV, and from the continuum region at a center-of-mass energy 
E cm = 10.52 GeV, just below the T(4S) resonance. The latter data set is used to subtract 
out the nonresonant continuum events produced at E cm = 9.46 GeV. The event sample taken 
at the T(1S) energy corresponds to an integrated luminosity of 62.5 pb _1 acquired during 
two different running periods. The sample of continuum events chosen for our background 
studies corresponds to an integrated luminosity of 91.3 pb _1 . 

To obtain a clean sample of hadronic events, we selected those events that had a minimum 
of three good charged tracks (to suppress contamination from QED events), a total visible 
energy greater than 15% of the total center-of-mass energy (to reduce contamination from 
two-photon events and beam-gas interactions), and an event vertex position consistent with 
the nominal e + e~ collision point to within ±5 cm along the e + e~ axis (z) and ±2 cm in 
the transverse (r — 0) plane. Backgrounds due to radiative Bhabha events with a converted 
photon (e + e~ — > e + e~7,7 — > e + e~) are reduced by requiring the total shower energy to be 
at least 15% of the available center-of-mass energy, but not more than 90% of the available 
center-of-mass energy, as well as by rejecting events that have thrust values approaching 1.0. 

Applying these cuts, we obtained 1.43 x 10 6 events from the two T(1S) data samples, 
collectively. 
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III. THE INCLUSIVE PHOTON SPECTRUM 



To obtain i? 7 , we first compiled an inclusive photon spectrum from the clusters of energy 
in the electromagnetic calorimeter. Only photons from the barrel region ( | cos ^ | < 0.7, 
where 9 1 is the polar angle of the shower) were considered. Photon candidates were required 
to be well separated from charged tracks and other photon candidates. The lateral shower 
shape was required to be consistent with that expected from a true photon. If the invariant 
mass of any two photon candidates fell within 15 MeV of the ir° mass, then both photons were 
rejected. Photons produced in the decay of a highly energetic 7r° would sometimes produce 
overlapping showers in the calorimeter, creating a so-called merged tt°. To remove this 
background, an effective invariant mass was determined from the energy distribution within 
a single electromagnetic shower. Showers whose effective invariant masses were consistent 
with those from merged 7r°'s were also rejected. Figure p] shows the inclusive spectrum that 
results from these cuts as a function of the scaled momentum variable, X 7 = p-y/E hea _ m . 

IV. BACKGROUND SOURCES 

The dominant source of background photons is asymmetric ir° decay. To remove this 
background, we developed a Monte Carlo simulation in which polar angle and event selec- 
tion effects were implicitly included. Modulo isospin breaking effects, one expects similar 
kinematic distributions between charged pions, which produce most of the charged tracks in 
T(1S) hadronic decays, and neutral pions, which produce most of the background pho- 
tons in the inclusive spectrum. By measuring the ratio of the true ir° momentum spectrum 
dN/dX^o to the charged track spectrum dN/dX n ±, the charged tracks themselves could then 
be used as a basis for simulating photons from 7r° decays. 

We therefore estimated the background due to photons produced in neutral meson decays 
as follows: for events that passed our selection criteria, we measured the ratio of efficiency- 
corrected 7r°'s to observed charged tracks as a function of momentum. Then, assuming 
that the angular distribution of 7r°'s is the same as that for charged tracks, the 3-momenta 
of the charged tracks were used to generate the expected background spectrum from n° 
decays (with the correct angular correlations implicit). The measured ratio provided the 
appropriate normalization. 

This approach had the advantage of being less model dependent than a Monte Carlo 
event generator, as the "generator" in this method was the data itself. It had the additional 
virtue that the absolute normalization of the 7r° background was simply determined by 
the number of accepted events. In addition to simulating the it — > 77 background, this 
technique was also used to account for r\ — > 77, to — > 7r°7, and rf — > j[p°, lu, 7] contributions. 
Figure [I| illustrates the corrected momentum spectra of these neutral mesons and charged 
tracks used to emulate their decays. 

Contributions from long lived neutral hadrons (neutrons, anti- neutrons, and K°'s) can 
also produce showers in the calorimeter. We used the LUND/JETSET 7.3 [|7| Monte Carlo 
simulation of T(1S) decays to estimate the number of long lived neutral hadrons in our event 
sample and a detector simulation based on the GEANT || package to determine how often 
these "residual showers" would pass the photon selection criteria. It was found that these 
hadrons represented a small contribution, not exceeding 3% for any value of X y . 
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A test of this background simulation method was performed using data collected from 
the continuum region, E cm = 10.52 GeV. Using a set of ratios for charged tracks to 7r°'s, 
77's, r) n s, and u's measured at this energy, we generated a photon spectrum and compared it 
to the inclusive spectrum from the continuum. With the exception of initial state radiation 
(whose contribution could be estimated from LUND and GEANT Monte Carlo simulations), 
the inclusive photon spectrum and simulated photon spectrum should agree. Figure |^ shows 
this comparison. We observe good agreement over the full range of A 7 .[] 

A. Subtractions and efficiencies 

This analysis has three major sources of background photons: (1) neutral hadrons (specif- 
ically, 7r°'s, 77's, r/"s, and cj's) produced in T(1S) decay, (2) neutral hadrons produced in non- 
resonant e + e~ — > qq'jX processes, and (3) radiative photons from the process e + e~ — > qq^y. 
By subtracting the dN/dX^ spectrum from the continuum data, scaled to correct for the 
differences in luminosity and cross section, we remove background from the latter two classes. 

The photon spectrum that we generated using charged particles collected at the T(1S) 
energy simulates the spectrum from the first two background classes combined, while the 
spectrum generated using charged particles from the continuum sample simulates only the 
second class. By subtracting these two spectra after appropriate scaling, we isolate the 
background spectrum of indirect photons from T(1S) decay. Hence subtracting the resulting 
spectrum from the data removes the first class of background. 

Figure f| shows the inclusive X 7 spectrum for data taken on the T(1S) resonance, with 
the different background contributions (non-resonant hadronic & radiative photons, resonant 
hadronic photons, and residual showers) overlaid. After subtracting these sources, what 
remained of the inclusive T(1S) spectrum was identified as the direct photon spectrum, 

T(1S)->7S0- 

To compare our data with predictions for the shape of the direct photon spectrum, we 
modified the theoretical distributions to account for attenuation and distortion from the 
finite detection efficiency and energy resolution. The most significant loss of direct photon 
events occurs at the high-X 7 region, arising from our requirement that an event have at least 
three good charged tracks. Unfortunately, hadronization in this kinematic regime is poorly 
understood. We considered two different T — > 'ygg Monte Carlo models with two different 
hadronization schemes and used a photon detection efficiency from the average of the two 
models (see Figure |5|). The difference in efficiency between the two models was incorporated 
into our systematic error. 

Trigger efficiencies have been evaluated directly from the data by determining the fraction 
of events passing a minimum-bias trigger. This efficiency, for all values of photon momentum 



Note: the Monte Carlo simulation of initial state radiation was not used as part of the final 
background subtraction. It is included in Figure || only to demonstrate that the background 
contribution to the inclusive photon spectrum is well modeled. Initial state radiation photons were 
automatically removed when we performed a scaled continuum subtraction to remove nonresonant 
contributions to the inclusive spectrum taken at E cm = 9.46 GeV. 
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considered in this analysis, exceeds 99%. 



V. COMPARISON WITH THEORETICAL MODELS 

A. Field Model 

As Figure ^ illustrates, the inclusive dN/dX 1 distribution increases rapidly in the low X 7 
region. This is due primarily to an overwhelming number of photons produced in tt° —>■ 77 
decays. However, to extract the total number of T(1S) — > ggj events and obtain i? 7 , we 
needed to integrate this spectrum along the entire scaled momentum axis. It was therefore 
necessary to rely on a model which fit well to that portion of the photon spectrum where the 
signal photons were clearly observable so that an extrapolation into the lower momentum, 
higher background region could be performed confidently. A number of attempts have been 
made to predict the shape of this spectrum @-|5| . In this analysis, we employed the model 
by Field || for our integration purposes. 

Figure || shows our photon spectrum with the background sources subtracted. To de- 
termine the number of direct photon events, N igg , from this spectrum, the data points in 
the region 0.30 < X 1 < 0.98 were fit to the modified (i.e., efficiency-attenuated and energy 
smeared) Field model; the only free parameter in this fit was the overall normalization. For 
comparison purposes, the modified lowest order QCD prediction,^ normalized to the same 
area as the Field model, has been overlaid. According to Field's model, about 85% of the 
direct photons that are produced in jgg decays lie within this portion of the momentum 
spectrum. According to our averaged detection efficiency of Figure ||, about 15% of those 
events are rejected by our shower and event-selection cuts. 

To determine the fraction of direct photons within our fiducial acceptance, we used a 
Monte Carlo simulation of the direct photon events, incorporating the QCD predictions of 
Koller and Walsh || for the photon angular distributions as a function of momentum. Ac- 
cording to their model, roughly 67% of the direct photons fall within our fiducial acceptance, 
\cos9\ < 0.7. Thus, our subtracted spectrum, within the limits of the fit and our fiducial 
acceptance, represents approximately 48% of the total direct photons produced in the T(1S) 
data sample. 

After integration of the fitted Field distribution in Figure |6| and corrections for finite 
acceptance, our data yields a total number of T(1S) — > •ygg decays, N igg = (2.652 ±0.038) x 
10 4 . 

To determine the number of three gluon events N ggg from the number of observed T(1S) 
hadronic events N^ S \ we first determined the number of continuum events under the T(1S) 
resonance (NjJ^P) from the observed number of T(4S) continuum events nJJ^P, accounting 
for the dependence of the cross-section on E^ m (= S): 

A£S S) = Njif .^m.^EL = 3.21 x 10 5 , (3) 

^cont <->T(lS) 



2 In the lowest order QCD prediction, the T system is treated as ortho-positronium decaying into 
three photons. 
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Next, we estimated the number of vacuum-polarization events N vp using the T(1S) — > 
branching fraction = 0.0248 fli"0" |, and Rr (is) — cr(e + e~ — > T(1S) — > qq)/a(e + e~ — 
Ti IS) 



/i 1 fi 



3.46 ±0.14 11 



N, 



vp 



-Rt(is) " 



N, 



T(1S) 
had 



1-35 



1.27 x 10- 



(4) 



From these values and Monte Carlo determined efficiencies for the various event types to 
pass our event selection selection criteria (see Table |I|), we determined 



N, 



mm 



T(1S) 
had 



-N. 



T(1S) 



cont 



• {j^-vp/ ^had) ^ygg ' (^7<7<?))A 



.'/.'/.'/ 



(9.657 ±0.010) x 10 5 . (5) 



From these values we obtained a value for R~: 



R^ 



N, 



imi 



N, 



2.75 ± 0.04% 



(6) 



.'/.'/.'/ 



B. Catani and Hautmann Modification to jgg Spectra 



Catani and Hautmann [ffl assert that in order to determine the total photon spectrum 
from T(1S) decays one must also consider fragmentation photons emitted from final-state 
light quarks produced in the initial heavy quarkonia decay. To properly measure a s , they 
claim, one must account for these additional photons, both in the shape of the spectrum, as 
well as in the QCD equations from which a s is extracted. They also provide a leading order 
estimate of the shape of the prompt photon spectrum due to this fragmentation component. 
In our analysis, we added this same component to the direct spectrum predicted by Field, 
modified the resulting spectrum for efficiency and energy resolution, and fit this distribution 
to our data using essentially the same method to determine N^ gg as described above. From 
that distribution, we measured i? 7 = (2.72 ± 0.04)% (see Figure |7|). Clearly, we do not 
yet have the requisite experimental sensitivity needed to verify the Catani and Hautmann 
model. 



VI. SYSTEMATIC ERRORS 

Table |I| summarizes the systematic errors studied in this analysis and their estimated 
effect upon i? 7 .[] The tracking efficiency and multiplicity modeling uncertainty was obtained 
by applying the two Monte Carlo models (with their different hadronization schemes) sep- 
arately, as opposed to their average. Including the 7r° veto reduces our statistical errors in 



3 Perhaps the largest potential errors arises from modeling the momentum and angular distribu- 
tions of the initial partons (i.e., the Field and Roller- Walsh predictions). Although these are not 
explicitly folded into our overall systematic error, it should be pointed out that our results are 
sensitive to these predictions. 
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the low X 1 region, but also adds to the uncertainty in our ability to accurately simulate 
this cut. The difference between the value of _R 7 obtained by applying the 7r° veto and the 
value obtained when we did not apply this veto constituted our second largest systematic 
uncertainty in i? 7 . By scaling the secondary photon spectrum by ±5%, we obtained the 
systematic error due to our uncertainty in the overall normalization of the secondary photon 
spectrum. We also compared results by using a different subtraction technique in which the 
non-resonant radiative contribution was subtracted using Monte Carlo simulated continuum 
events generated at the T(1S) center-of-mass energy. This allowed us to extract a value 
of i? 7 independent of any non-resonant data taken at energies other than 9.46 GeV. The 
estimated uncertainty in the number of three gluon events can be directly translated to an 
uncertainty in i? 7 . To check against possible systematic effects due to different running 
conditions, we analyzed the two T(1S) data samples separately. Finally, we included the 
total error (statistical and systematic, combined in quadrature) quoted by ARGUS in their 
measurement of the ratio of hadronic to muonic cross-sections, Rytis) ■ 

Table |II| compares the results of this analysis with those obtained by previous exper- 
iments in which the observed number of T(1S) — > jgg events were also determined using 
Field's model. 



VII. EXTRACTION OF QCD PARAMETERS 



We now relate the value of _R 7 to the fundamental QCD parameters which we wish to 



measure, following Sanghera | 15fl . 

The decay width T — > jgg has been calculated by Lepage and Mackenzie in terms 
of the coupling strength at the energy scale characterizing this decay process, a s (Mx): 



r(T- T 00) 8(tt 2 -9) 



a 2 a {M- 



1 + (3.7 ±0.4) 



«s(M T ) 



7T 



r(T -> fx+fi ) 9na QED 
Expressing this ratio in terms of a leading-order power series in a s (/z), we have: 



(7) 



7T 



7T 



l-nb^ln 



+ (3.7 ±0.4) 



(8) 



where A, = b = (33-2n / )/12 7 r, and n/ is the number of light quark flavors which 

participate in the process (n/ = 4 for T(1S) decays). 



Similarly, the decay width T — > ggg has been calculated by Bardeen et al. |17j and 
expressed by Lepage et al. 0,11^] as: 



T(T^ggg) _ 10(vr 2 - 9) a s 3 (M x ) 



T(T -> fi+fi- 



8lnel 



QED 



a s {M r ) 



[2.770(7)/5 - 14.0(5)] + • 



(9) 



where (3q = 11 — (|)^/, and e b = — |, the charge of the b quark. Again, we can express this 
in terms of the renormalization scale: 



T(T -> ggg) 

r(T - 



A r 



a s (n) 



7T 



7T 



3nboln 



(JL. 



-)B f n f + B i 



(10) 
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where A Q = 10 ^V 9) / , B f = 2.770 ± 0.007, and B t = 16.47 ± 0.58. 

a 8ie (, a QED 

The strong coupling constant a s can be written as a function of the basic QCD parameter 
A.jjg, defined in the modified minimal subtraction scheme [10| , 



„(„) = ! ( 1 h 'ndn^/A^m 



153-19^ 
24?r 2 



where b\ 

Note that the scale dependent QCD equations (§) and (pTOf ) are of finite order in a s . If 
these equations were solved to all orders, then they could in principle be used to determine 
Rry independent of the renormalization scale. Because we are dealing with calculations that 
are of finite order, the question of an appropriate scale must be addressed. 

The renormalization scale may be defined in terms of the center of mass energy of the 
process, fi 2 = f^E 2 m , where f ^ is some positive value. But QCD does not tell us a priori 
what fn should be. One possibility would be to define \i = E cm ; that is / M = 1. A 
number of phenomenological prescriptions |1^JT8|J20|J21[| have been proposed in an attempt 



to "optimize" the scale. However, each of these prescriptions yields scale values which, in 



general, vary greatly with the experimental quantity being measured fl5| . 

In this analysis, we have determined Aj^g over a range of scale values. This was done by 
comparing our measured value of _R 7 with the ratio of equations (El) and ( |i~0"|) in which a s 
was replaced by the expression in equation (|ilf ), thereby providing a relationship between 
Ajjg, fp, and i? 7 . Thus for each assumed value of / M , Ajjg was numerically determined 
as a function of i? 7 . The resulting Ajjg versus / M dependence is shown in Figure ||. This 
dependence was parameterized by the form 

^Ws(f) = AmK/o) + ( c i)^(^) + ( c i + C2 ^ 2 (|) ^ 12 ) 

where fo is the value of around which Aj^ is minimally dependent on / M , given by 
(dAjjg/df^ = 0). In this analysis, we determined fo = 0.107, Aj^g(f ) = 168.62 MeV, 
d = 7.74, c 2 = 14.68. 

By parameterizing the results of the analysis in this manner, one can easily extract QCD 
parameters at any scale within the range of the parameterization, 0.10 < < 1.0, and 
compare with other results. For example, the mean value between Ajjg(f^ = 0.107) (where 
Ajys is a minimum), and A^gg(/ M = 1.0), is 233 ± 11 ± 59 MeV. The uncertainty of the 
parameterization due to the theoretical uncertainties of the parameters in equations (||) 
and fllCf ) has been included in the systematic error of A^r. Substituting this value for A-jg-g 
into equation (|Tl~D, and using /i = Mms) we find for a s : 

tts(M T(ls) ) = 0.163 ± 0.002 ± 0.009 ± 0.010, (13) 

where the additional error of ±0.010 arises from the difference, about 64 MeV, between the 
mean value of Ajjs an d the values at each of the parameterization limits, / M = 0.107 and 
fn = 1.0. Extrapolating this result to fi = Mz, and assuming continuity of a s across the 
5 flavor continuum threshold |22| (which implies that A-^- is a step function across the 5 
flavor threshold), we obtain from equation (|TTD: 
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a s (M z ) = 0.110 ± 0.001 ± 0.004 ± 0.005. 



(14) 



This result is lower, although in acceptable agreement with the average value of a s (Mz) = 
0.122 ± 0.007 presently quoted by the Particle Data Group [|Kj]. It is worth noting that the 
value of Ajj-g obtained by previous experiments Hl3j|l4ll using the fixed-scale procedure of 



ref. |18[ is in agreement with our value at fo where the scale dependence on Ajrq becomes 



minimal. 
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TABLE I. Monte Carlo Determined Event Efficiencies 



Event Type Symbol Efficiency 

three gluon e ggg 0.9938 

'generic' T(1S) hadronic thad 0.9985 

vacuum polarization e vp 0.9480 

direct photon e 799 0.9419 



TABLE II. Systematic Errors 

Uncertainty Source 6FLy(%) 

Tracking efficiency and multiplicity modeling 0.12 

7T° veto 0.07 

continuum subtraction 0.04 

Eggg 0.03 

pseudo-photon spectrum 0.03 

Luminosity and E cm scaling 0.02 

T(1S) data samples used separately 0.02 

<Hfcr(is) 0.01 



TABLE III. Comparison with Other Experiments 

Experiment i? 7 (%) 

CLEO 1.5 @ 2.54 ±0.18 ±0.14 

ARGUS || 3.00 ± 0.13 ± 0.18 

Crystal Ball |L§ 2.7 ± 0.2 ± 0.4 

This measurement 2.75 ±0.04 ±0.15 
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X 7 

FIG. 1. The inclusive distribution of candidate photons as a function of scaled momentum 
Xy = p^/Ebeam, from data taken at the T(1S) center-of-mass energy. Note the log scale of the 
y-axis. 
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FIG. 2. Efficiency corrected 7r° momentum spectrum, Monte Carlo generated 77, r/, and u; spec- 
tra, and observed charged tracks' momentum spectrum as a function of scaled particle momentum, 
X w = ^/-^beam- I n this notation, V refers to any of the neutral spectra or charged tracks. 
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FIG. 3. A comparison of the inclusive photon spectrum from continuum data taken below the 
T(4S) resonance (histogram) with the simulated background spectrum from 7r°'s, r/'s, r/'s, and w's 
produced by non-resonant processes at this energy. To illustrate the magnitude of the initial state 
radiative (ISR) correction, the simulated spectrum both with (dark squares) and without (open 
circles) the (Monte Carlo determined) radiative contribution are overlaid. 
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FIG. 4. The inclusive X 7 spectrum (histogram) for data taken at the T(1S) resonance, along 
with background contributions due to non-resonant processes (open circles), resonant hadronic 
decays (dark squares), and other residual sources (inverted triangles). 
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FIG. 5. The direct photon detection efficiency as a function of scaled photon momentum X 
determined by averaging the two Monte Carlo simulation models (see text). 



19 



1500 



3031096-009 



t — i — r 



1000 



500 - 



i — i — i — r 



t — i — i — r 



i — i — i — r 



Background Subtracted Data 
Modified Field Model 
(Fitted Region) 
Modified Field Model 
(Extrapolated Region) 

Lowest Order QCD Prediction 




0.25 



0.50 

X 7 



J I I L 

0.75 1.00 



FIG. 6. The background subtracted (i.e., direct) photon spectrum (dark squares). The data 
points in the region 0.3 < X-y < 0.98 are fit to Field's model (histogram). The only free parameter in 
the fit was the overall normalization. The fit to Field's model, extrapolated into the low momentum 
region, is also shown (dashed line) as well as the lowest order QCD prediction (dotted line) over the 
full kinematic regime, which has been normalized to the same area as the Field model. The errors 
shown are purely statistical. Data points in the region X 7 < 0.3 appear systematically shifted 
above the Field line, however, a one sigma shift in the magnitude of the simulated background line 
in Figure |I| would drastically alter the distribution at the low end. Systematic errors, had they 
been included on these points, would be well off scale. 
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FIG. 7. The background subtracted photon spectrum (dark squares) fit to the Field distribution 
with the added fragmentation component predicted by Catani and Hautmann. Again, the errors 
shown on the data points are purely statistical. 
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FIG. 8. Ajgg as a function of scale parameter f^, as governed by the functional dependence of 
i? 7 on Aj^ and / M . 
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